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Abstract. 

We discuss shell-model calculations based on the use of low-momentum interactions derived 
from the free-space nucleon-nucleon potential. A main feature of this approach is the 
construction of a smooth potential, Viow-k, defined within a given momentum cutoff. As a 
practical application of the theoretical framework, we present some selected results of our current 
study of nuclei around doubly magic ^"^^Sn and ^"^^Pb which have been obtained starting from 
the CD-Bonn potential. Focusing attention on the similarity between the spectroscopy of these 
two regions, we show that it emerges quite naturally from our effective interactions without use 
of any adjustable parameter. 



1. Introduction 

In the last decade, shell-model calculations employing realistic effective interactions derived from 
modern nucleon-nucleon (NN) potentials have entered the main stream of nuclear structure 
theory [T]. As is well known, the first problem one is confronted with in this kind of calculations 
is the strong short-range repulsion contained in the bare A^A^ potential Vat at, which prevents its 
direct use in the derivation of the shell- model effective interaction Vqs. The most popular way 
to overcome this difficulty has long been the Brueckner G-matrix method. However, a few years 
ago a new approach [2j was proposed which consists in deriving from Vnn a low-momentum 
potential, Vjow-k) that preserves the deuteron binding energy and scattering phase shifts of Vnn 
up to a certain cutoff momentum A. This is a smooth potential which can be used directly to 
derive 14fj, and it has been shown [21 13] that it provides an advantageous alternative to the use 
of the G matrix. In this connection, it should be mentioned that Vjow-k potentials are currently 
being used in various nuclear theory contexts, such as the study of few-body systems and no-core 
shell-model calculations [U E] . 

Making use of the Mow-k approach, we have recently studied [U [3 El [9] several nuclei beyond 
doubly magic ^^'^Sn, showing that their properties are well accounted for by a unique shell- 
model Hamiltonian with single-particle energies taken from experiment and two-body effective 
interaction derived from the CD-Bonn A^A^ potential [10]. 

Motivated by the very good results obtained in the ^^^Sn region and by the existence of a 
close resemblance [HI [I2l [I3l [H] between the spectroscopy of this region and that of nuclei 
around stable ^''^Pb, we have found it challenging to perform a comparative study of these two 



regions [E]- In this paper we present some results from this study, focusing attention on the 
proton-proton, neutron-neutron and proton- neutron multiplets in the three far-from-stabihty 
nuclei ^^"^Te, ^^^Sn, ^^^Sb and in their counterparts in the ^°^Pb region, ^^°Po, ^^°Pb, and ^^°Bi. 

We start by giving an outline of the theoretical framework in which our shell-model 
calculations are performed and then present and discuss our results. A short summary is given 
in the last section. 



2. Theoretical framework 

In the framework of the shell model an auxiliary one-body potential U is introduced in order to 
break up the nuclear Hamiltonian, written as the sum of the kinetic term and the A^A^ potential, 
into a one-body component Hq, which describes the independent motion of the nucleons, and a 
residual interaction Hi. Namely, 

A 2 

H = Y.^ + T.^ij = T + Vnn = {T + U) + {Vnn -U) = Ho + Hi. (1) 

1=1 i<j 

A reduced model space is then defined in terms of the eigenvectors of Hq and the 
diagonalization of the original Hamiltonian H in an infinite Hilbert space is reduced to the 
solution of an eigenvalue problem for an effective Hamiltonian i^eff in a finite space. 

The Hamiltonian i^eflf can be derived by way of the Q-box folded-diagram expansion (see 
Ref. [T]). This implies as first step the calculation of the so-called Q-box, which is made up 
of an infinite collection of irreducible and valence-linked Goldstone diagrams. Once the Q- 
box has been calculated at a given order, the infinite series of the folded diagrams has to be 
summed up. From this procedure an effective Hamiltonian is obtained containing both one- and 
two-body components. Usually, only the two-body term V^g is retained, while the one-body 
contributions, representing the theoretical single-particle energies, are subtracted and replaced 
with single-particle energies taken from experiment [16j . 

However, as mentioned in the Introduction no modern NN potential can be used in a 
perturbative nuclear structure calculation, unless its strong repulsive core is firstly "smoothed 
out" . Here, we do not embark on any discussion of how the effective Hamiltonian is derived, and 
refer to [1], where a detailed description of the whole procedure can be found. Rather, in the 
following we focus on the Mow-k approach. We first outline the essential steps for the derivation 
of the Mow-k based on the Lee-Suzuki similarity transformation method [TT], and then describe 
its main features. 

Let us consider the similarity transformation on the Hamiltonian ^ 

n = X-^HX, (2) 

where the operator X is defined in the whole Hilbert space. We now introduce a cutoff 
momentum A that separates fast and slow modes to the end of deriving from the original 
V/Viv a low- momentum potential satisfying a decoupling condition between the low- and high- 
momentum spaces. 

The low-momentum space is specified by 

P = J dp\p){p \ , p<A (3) 
where p is the two-nucleon relative momentum, and the decoupling equation reads 

QHP = 0, (4) 

with Q = 1 — P being the complementary fast-mode space. The low-momentum Hamiltonian is 
then given by 



i7io„-k = PHP, (5) 

and it can be easily proved that its eigenvalues are a subset of the eigenvalues of the original 
Hamiltonian. 

There are, of course, different choices for the transformation operator X. We take 

X = e-, (6) 
where the wave operator oj satisfies the conditions: 

uj = QujP, (7) 
PujP = QojQ = PujQ = 0, (8) 

the former implying that 

X = l + u;. (9) 
From Eq. ([5]) the low-momentum potential Viow-k can be defined as 

Mow-k = ^low-k - PTP. (10) 

Employing transformation ([9]), this equation is written as 

Mow-k = PVnnP + PVnnQuj, (11) 

while Eq. (jj]) becomes 

QVnnP + QHQiu - LoPHP - loPVnnQoo = 0. (12) 

The solution of Eq. (fO]) gives the value of uj needed to obtain Fiow-k- 

This decoupling equation can be solved by means of the iterative technique for non-degenerate 
model spaces proposed in ^H], which is now sketched. We define the operators: 

p{uj) = PHP + PHQlo, (13) 
qiio) = QHQ-uPHQ, (14) 

in terms of which one can write 

xo = -{QHQ)-^QHP , 
xi = g(xo)~"^xop(xo) , 

Xn = q{xo-\-Xi + ...+Xn-l)~^Xn-lp{xo + Xi + ... + Xn-l) ■ (15) 

Once the iterative procedure has converged, Xn 0, the operator lj is given by 

n 

^n = ^Xi, (16) 
i=0 

In applying this method, we have employed a momentum-space discretization procedure making 
use of an adequate number of Gaussian mesh points [19] . 



It is worth mentioning that the equation for Vjow-k obtained using the similarity 
transformation of Lee and Suzuki is the same as that one can derive from the T-matrix 
equivalence approach [H [2] . This means that the obtained low-momentum potential not only 
preserves the deuteron binding energy given by the original A'^A^ potential, but also its low- 
momentum (< A) half-on-shell T matrix. 

The above Viow-k is however not Hermitian, which is not convenient for various applications, 
as for instance its use in the derivation of shell- model effective interactions. This V|ow-k may 
be transformed by means of the familiar Schmidt orthogonalization procedure, which leads to 
a Hermitian new T^o„_k. As suggested in |18j . another transformation, based on the Cholesky 
decomposition of a symmetric and positive definite matrix, can be used to this end. In fact, the 
matrix P(l -\-lo~^uj)P, being symmetric and positive definite, admits this decomposition, 

P(l + uj+uj)P = PLL^P, (17) 

where L is a lower triangular matrix and its transpose. Since L is real matrix defined within 
the P-space, we may write our transformation as 

Z = L^, (18) 

and the corresponding Hermitian Vjow-k is 

V^t°ik = PL'^Pi.Ho + Vio^_y,)P{L-YP - PHoP (19) 

In Refs. |20l [Tj, it has been shown that this Hermitian interaction, as well all the family of 
Hermitian interactions which can derived from Mow-k by means of different transformations, 
preserve the full-on-shell T matrix, and consequently the phase shifts of the original Vnn- 

The so-obtained V|ow-k is a smooth potential that can be used directly within the Q-box 
folded diagram theory to derive the shell-model effective interaction. Actually, it represents an 
advantageous alternative to the G-matrix approach owing to the fact that it does not depend 
either on the energy or on the model space. This is at variance with the G matrix, which is 
defined in the nuclear medium. In this connection, it is worth mentioning that the merit of 
the Mow-k within the context of realistic shell-model calculations has been assessed by several 
studies evidencing that Mow-k results are as good, or even slightly better than, the G-matrix 
ones [21 [3]. 

Finally, it is a remarkable feature of the Mow-k approach that different A^A^ potentials lead 
to low-momentum potentials, which are quite similar to each other [H 121] . 

3. Two-valence particle nuclei in ^^^Sn and 208p]^ regions 
3.1. Outline of calculations 

In this paper, we present some results of our current shell-model study of nuclei with two valence 
nucleons in the ^^^Sn and 208p]-, regions, which have been obtained starting from the CD-Bonn 
potential renormalized by use of a V|ow-k with a cutoff momentum of A = 2.2 fm~^. 

In our calculations for ^^^Te, ^^^Sn, and ^^^Sb we assume that the the valence protons occupy 
the five levels 0(77/2, lc^5/2) 1^^3/2, "^^1/2^ and 0/iii/2 of the 50 — 82 shell, while for the neutrons 
the model space includes the six levels 0/^9/2, I/7/21 I/5/21 2^3/2, 2^1/2, and 0^13/2 of the 82 — 126 
shell. Similarly, for ^^'^Po, ^^"Pb, and ■^^'^Bi we take as model space for the valence protons 
the six levels of the 82 — 126 shell and let the valence neutrons occupy the seven levels 159/2; 
0^11/2; 0ii5/25 2(^5/2, 3si/2, 1^7/2 and 2^3/2 of the 126 — 184 shell. As regards the adopted single- 
particle neutron and proton energies, they can be found in Refs. [6] and [22] for ^^^Sn and ^ospb, 
respectively. 
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Figure 1. (a) Proton-proton (1^/7/2)'^ multiplet in ^^^Te.(b) Proton-proton (2^59/2)^ multiplet 
in ^^'^Po. The theoretical results are represented by open circles while the experimental data by 
solid triangles. 

As mentioned in section 2, the two-body matrix elements of the effective interaction are 
derived within the framework of the Q-box folded-diagram expansion. We include in the Q- 
box all diagrams up to second order in the interaction, given by the Mow-k potential plus the 
Coulomb force for protons. These diagrams are computed within the harmonic-oscillator basis 
using intermediate states composed of all possible hole states and particle states restricted to 
the five proton and neutron shells above the Fermi surface. The oscillator parameter is 7.88 
MeV for A = 132 region and 6.88 MeV for the A = 208 region, as obtained from the expression 
huj = 45 A~^/^ — 25 A^"^/^. The calculations have been performed by using the NUSHELLX 
code 123]. 

3.2. Results 

In figures 1, 2, and 3, we present the experimental [MIES] and calculated excitation energies of 
the lowest states in the nuclei with two-proton, two-neutron, and one proton-one neutron beyond 
doubly magic ^^^Sn and ^''^Pb. We see that the agreement between theory and experiment is 
very good for all six nuclei considered, the discrepancies being well below 100 keV for most of 
the states. 

All the calculated states reported in these figures are dominated by a single configuration, 
whose percentage ranges from 80% to 100%. In particular, they correspond to the proton-proton 
multiplets {1^97/2)'^ and (vr/i9/2)^ in ^^^Te and ^-'^'^Po, to the neutron-neutron multiplets (1^/7/2)^ 
and (j^fi'9/2)^ in ^^'^Sn and ^-'^'^Pb, and to the proton-neutron multiplets t^Qt h j2 and vr/i9/2J^59/2 
in i34sb and ^lOfii. 

These figures evidence the striking resemblance between the behavior of the multiplets in the 
three pairs of counterpart nuclei, ^^^Te and ^^°Po, ^^'^Sn and ^^''Pb, ^^^Sb and ^^°Bi. We may 
only note that the curves relative to the ^'^^Pb neighbors are generally located slightly below 
those for the counterpart nuclei in ^^^Sn region. In particular, from figures 3a and 3b we see 
that the two proton-neutron multiplets show a sizable energy gap between the 2~ state and the 
nearly degenerate 0~ and 1~ states, as well as a distinctive staggering, with the same magnitude 
and phase, between the odd and the even members starting from the 3~ state. As regards the 
two-identical-particle multiplets, figures 1 and 2 show four curves having all the same shape. 
It is worth noting, however, that the curves for the two-valence-proton nuclei are located in an 
energy interval larger than that pertaining to the two-valence-neutron nuclei. This means that 
in both ^^^Sn and ^osp]^ regions a weakening of the pairing gap exists for nuclei with two- valence 
neutrons with respect to those with two- valence protons. 




Figure 2. (a) Neutron-neutron ij^g'jj'})'^ multiplet in ^^^Sn.(b) Neutron-neutron (vr/ig/2)^ 
multiplet in ^^'^Pb. The theoretical results are represented by open circles while the experimental 
data by solid triangles. 



It is worth mentioning that the resemblance between Sn and Pb regions was first 
pointed out by Blomqvist[TT], who noticed that every ^^^Sn single-proton and -neutron level, 
characterized by quantum numbers (n/j), has its counterpart around ^'^'^Pb with quantum 
numbers {nl + Ij + 1). However, until recent years the scarcity of information for nuclei around 
^^^Sn, which lies well away from the stability line, has prevented a detailed comparative study 
of the two regions. Nowadays, new data have become available which support the similarity 
between their spectroscopies. In our calculation, this similarity emerges quite naturally from 
our effective interaction which we have derived from a realistic NN potential. 
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Figure 3. (a) Proton-neutron vr (77/2^/7/2 niultiplet in ^^^Sb.(b) Proton-neutron vr/i9/2 1^(79/2 
multiplet in ^^'^Bi. The theoretical results are represented by open circles while the experimental 
data by solid triangles. 



4. Summary 

We have briefly discussed here the theoretical framework for realistic shell-model calculations 
wherein use is made of low-momentum interactions derived from the free NN potential. We have 
shown how a smooth low-momentum potential V|ow-k can be constructed, which preserves the 
deuteron binding energy and scattering phase shifts of the original Va^at up to a given momentum 
cutoff. We have then presented the results of a shell-model study of nuclei around doubly magic 



Sn and Pb, focusing attention on proton-proton, neutron-neutron and proton-neutron 
multiplets. The results obtained for the three nuclei ^^^Te, ^^'^Sn and ^^^Sb have been compared 
with those for ^^'^Po, ^^'^Pb and ^^'^Bi, which are their counterparts in the region of ^''^Pb. In 
both low-momentum effective interaction derived from the CD-Bonn A^A^ potential has 

been employed. It should be stressed that no adjustable parameter appears in our calculations. 

Our results for all six nuclei are in very good agreement with the experimental data and 
account for the striking resemblance between the behavior of the multiplets in the ^^^Sn and 
208p]2 regions. This stimulates further studies to find out whether this resemblance extends 
beyond the two-valence-particle nuclei. 
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